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ABSTRACT 

To obtain a quantative understanding of the cratering process over a broad range of 
conditions, we have numerically computed the evolution of impact induced flow fields and 
calculated the time histories of the major measures of crater geometry (e.g. depth diameter, 
lip height. . . ) for variations in planetary gravity (0 to 10 9 cm/s 2 ), material strength (0 to 
140 kbar), thermodynamic properties, and impactor radius (0.05 to 5000 km). These 
results were fit into the framework of the scaling relations of Holsapple and Schmidt 
[1987]. 

We describe the impact process in terms of four regimes: (1) penetration, (2 ) 
inertial, (3) terminal and (4) relaxation. During the penetration regime, the depth of 
impactor penetration grows linearly for dimensionless times x = (Ut/a) <5.1. Here, U is 
projectile velocity, t is time, and “a” is projectile radius. In the inertial regime, X >5.1, 
the crater grows at a slower rate until it is arrested either by strength or gravitational forces. 
The crater depth, D, and diameter, d, normalized by projectile radius is given by 
D/a =1.3 (Ut/a) 0 - 36 and d/a= 2.0= (Ut/a) 0 - 36 . For strength dominated craters, growth 
stops at the end of the inertial regime, which occurs at (Yeff/pU 2 )-^- 78 , where Y e ff is the 
effective planetary crustal strength. In gravity dominated craters, growth stops when the 
gravitational forces dominate over the inertial, which occurs at x = 0.92 (ga/U 2 )* 0 - 61 . In 
strength dominated craters, the growth stops at the end of the inertial regime, which occurs 
at x = 0.33 (Yeff/pU 2 )* 0 - 78 . In the gravity and strength regimes the maximum depth of 
penetration is Dp/a =1.2 (ga/U 2 )- 0 - 22 and Dp/a = 0.84 (Y/p U 2 ) -0 - 28 , respectively. 

The transition form simple bowl-shaped craters to complex-shaped crater is found 
to result from the transition from strength dominated to gravity dominated craters. The 
diameter for this transition to occur is given by dj = Y/pg, and thus scales as g* 1 for 
planetary surfaces where strength is not strain-rate dependent. This result agrees with 
crater-shape data for the terrestrial planets [Chapman and McKinnon, 1986]. 
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We have related some of the calculable, but non-observable parameters which 
characterize the cratering process ( e.g. maximum depth of penetration, depth of 
excavation, and maximum crater lip height) to the crater diameter. For example, the 
maximum depth of penetration relative to the maximum crater diameter is 0.58, for strength 
dominated craters, and 0.28 for gravity dominated craters. These values imply that 
impactors associated with the large basin impacts penetrated relatively deeply into the 
planet's surface. This significantly contrasts to earlier hypotheses in which it had been 
erroneously inferred from structural data that the relative transient crater penetration 
decreased with increasing diameter. Similarly, the ratio of the maximum depth of 
excavation relative to the crater diameter is a constant s 0.1, for gravity dominated craters, 
and ^).2 for strength dominated craters. This result implies that for threshold impact 
velocities less than 25 km/sec, where significant vaporization begins to take place, the 
excavated material comes from a maximum depth of 0.1 times the transient crater diameter 
Moreover, we find apparent final crater diameter is approximately twice the transient crater 
diameter at diameters of greater than dp8.6 Y/pg. 
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INTRODUCTION 

The recent decades of planetary observation and exploration have lead to the 
conclusion that the impact of solid bodies is one of the fundamental processes in the origin 
and evolution of the solar system.While the impact process is conceptually easy to 
visualize, the detailed quantitative description of the mechanics has been difficult and 
illusive (e.g. Melosh [1989]). The overall objective of this study is to obtain a quantative 
understanding of the planetary cratering process over a broad range of conditions. 
Specifically the objectives are to establish quantitative scaling relationships for (1) the 
temporal evolution of the key measures of the crater geometry (depth, diameter, and lip 
height), (2) the maximum depth of penetration and excavation, and (3) the transition 
from simple to complex craters. We address these objectives over a broad range of 
planetary gravities, material strengths, and impactor sizes (from meters to those that formed 
the multiringed basins). 

APPROACH AND SCOPE 

We modeled the normal impact of spherical projectiles on a semi-infinite planet 
surface over a broad range of conditions using numerical techniques. We do not address 
the effect of a planet’s atmosphere upon the cratering process in the present paper and refer 
readers to [Melosh, 1989 Chapt. XI; O'Keefe and Ahrens, 1988; Roddy et al., 1987; 
Schultz and Gault, 1981]. To calculate the impact-induced flow fields within the solid 
planets, we used the Eulerian-Lagrangian code developed by Thompson [1979]. The key 
equation of state parameters along with the mechanical parameters that were used in the 
code to represent a typical silicate impactor and planet are listed in Table 1. The constitutive 
model accounted for elastic-plastic hydrodynamic responses [Thompson, 1979]. Other 
models are being examined and will be reported upon in the future. We varied some of the 
equation of state parameters such as the melt enthalpy and vaporization enthalpy, and also 
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the mechanical properties such as the yield strength. In all cases, the material properties of 
the impactor and planet were identical. 

The impact parameters and variables are listed in Table 1 along with the range over 
which they were varied. We nondimensionalized these using the formalism of Holsapple 
and Schmidt [1987]. The magnitudes of the four dimensionless parameters ga/U 2 , 
Y/pU 2 , Hm/U 2 and H va pAJ 2 are measures of the dominant mechanisms controlling the 
cratering process. We varied these parameters over a range of conditions so as to 
determine when each of these was the dominant parameter that described the impact 
process. 

The inverse Froude number, ga/U 2 , is a measure of the gravitational forces relative 
to the pressure forces. The inverse Cauchy number, Y/pU 2 , is a measure of the planetary 
strength relative to the pressure forces; and the melt, Hm/U 2 , and vaporization, HvapAJ 2 , 
numbers are measures of the relative importance of melting and vaporization. In 
determining the range of the inverse Froude number, we had a choice of varying either the 
projectile radius or the gravitational acceleration; for ease of computation, we fixed the 
radius and varied the gravitational acceleration over six orders of magnitude. This is 
equivalent to varying the impactor radius from 5 m to 5000 km. In determining the range 
of the vaporization number, HyapAJ 2 , we restricted the impact velocities to 12 km/s so as to 
not get into the impact regime where there are significant vaporization effects. This limits 
the validity of the scaling laws to velocities less than 30 km/s for typical silicates. We have 
studied the effects of high speed impacts (>30 km/s) and will report in detail on these 
elsewhere [O’Keefe and Ahrens, 1989]. 

The results of the computations represent a very large amount of computer output 
with over -10 4 variables being calculated for each time step. In this paper we will report 
on the geometrical measures of the cratering process such as the depth, diameter and lip 
height The depth is defined as the distance from the initial planetary surface and the 
planetary surface at a given time at the centerline of the impact; the diameter is defined as 
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the distance between the interface between the impactor and planetary surface at the initial 
planetary surface height; the lip height is measured relative to the original planetary surface 
level (e. g. see figure 1). The dimensionless depth and diameter histories are summarized 
in Table 2. In addition, we have included a series of detailed flow field plots; these 
represent the development of a simple bowl shaped crater (figure 1) and a complex crater 
exhibiting central peak and ring formation (figure 2) .To determine the displacement of 
planetary material during the cratering process, we placed massless tracer particles at 
various positions and computed their trajectories which are also shown. 

CRATER SCALING REGIMES 

The cratering process can be described in terms of at least four temporal regimes; 
penetration, inertial, terminal and relaxation. Schmidt and Holsapple called these first three 
regimes the early stage, intermediate and late. The numerical approach taken here can be 
used to describe the first three; the fourth, the terminal regime is the result of long term 
equilibration and requires a different analytical approach such as viscous relaxation 
methods (e.g. Melosh [1989], Chapt. 8). In the following section we will describe the first 
three regimes and determine the numerical values of the parameters in the Holsapple and 
Schmidt [1987] scaling laws. We will describe the depth, diameter, and other geometrical 
measures in that order. The scaling laws and the numerical values of the parameters derived 
from the present calculations are listed in Table 3. 

PENETRATION REGIME 

The penetration regime is characterized by the transfer of the kinetic energy of the 
projectile to the planetary surface. In the case where the impactor and planetary surface 
have similar properties, about half of the impactor kinetic energy is transferred to the planet 
during the penetration regime [O'Keefe and Ahrens, 1977]. During this time, the impactor 
drives a strong shock wave into the planet and also deforms and lines the growing crater 
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cavity wall (Figures la,b). The material properties of the impactor and the planet and the 
planet’s gravitational acceleration are not important during this time, except for relative 
differences in density between the impactor and the planet [Holsapple and Schmidt, 1987]. 
Depth 

The crater depth grows linearly with time in the penetration regime. Shown in 
Figure 3 is a plot of the evolution of the interface between the impactor and the planet at the 
centerline of impact for a wide range of impact conditions. Note that all the cases fall on a 
single line independent of strength or gravitational acceleration. This linear growth is 
expected during the early times when the shock at the interface is nearly planar and is 
predicted from planar shock wave theory. The depth(D)/impactor radius(a) evolves as 

D/a = j (Ut/a) (1) 

where U is the impact velocity. The factor of two comes from the fact that upon impact of 
bodies of like impedances, the interface velocity is one-half the impactor velocity. 

The depth grows linearly with time until it enters the inertial regime This was 
predicted by Holsapple and Schmidt [1987] to occur at dimensionless time, x, of 5.1 for 
like impacts and this is substantiated in Figure 3. 

Diameter 

The diameter does not have a simple growth law during the penetration regime. The 
impact of a near spherical impactor will produce small amounts of jetting and vaporization 
near the initial point of impact even at low velocities; the magnitude depends sensitively on 
the details of the impactor-planetary surface geometry [KiefFer, 1976]. The penetration 
regime for diameter ends when the maximum lateral extent of the impactor penetrates the 
surface which, from geometrical considerations for a spherical impactor, occurs at 
dimensionless time of 

x=Ut/a = 2 (2) 
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Shags 


The shape of the transient cavity changes rapidly during the penetration regime. The 
depth grows more rapidly than the diameter during the penetration time for the depth, and 
the scaling of the evolution of the diameter to depth is given from the ratio of equations 1 
and 6 

d/D = 4.0 ( Ut/a yM* ; 1 < ( Ut/a ) < 5.1 (3) 

INERTIAL REGIME 

The inertial regime is characterized by the quasi-hemispherical expansion of the 
crater cavity. In this regime, the geometry of the cavity does not change, and the projectile 
is deformed into a thin hemispherical shell that lines the transient crater cavity. The strong 
shock that was attached to the projectile during the penetration time has now propagated 
away from the the cavity region [Bjork, 1961; O'Keefe and Ahrens, 1977]. Bjork called 
this regime the detached shock regime. As in the penetration regime, the material strength 
properties and gravity do not play a role in the evolution, however, the termination of this 
regime occurs when either strength or gravitational effects arrest crater growth. 

Depth 

The depth grows as a simple power law in the inertial regime and it is independent 
of either strength or gravitational acceleration (see Figure 3). A similar result for the depth 
was found by Holsapple and Schmidt [1987 Figure 4] for a variety of impact conditions. 
Eq.4 is obtained by fitting the depth versus time curve (Fig. 3) in the inertial regime. The 
scaling law is given by 

D/a = 1.3 (Ut/a) 0 - 36 (4) 
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The magnitude of the inertial regime exponent (S=0.36) is related to the coupling 
coefficient of Holsapple and Schmidt [1987] by the relation 

\l = S/(l-S) (5) 

Eq. 5 gives p = 0.56, which is consistent with the range of values found by Holsapple 
and Schmidt [1987] for a range of impact conditions. In the following sections all of the 
scaling laws will have exponents that are functions of p and thus are related to the exponent 
in the inertial regime. 

Di am eter 

The diameter also grows as a simple power law in the inertial regime (e.g. Figure 
4). The scaling for the diameter from a fit to the results in Fig. 4 is given by 

d/a = 2.0 (Ut/a) 0 - 36 (6) 

Note that the exponent is the same for both the depth and diameter. 

Shape 

The shape is similar throughout the duration of the inertial regime. The ratio of the 
diameter to depth is 

d/D = 1.6 : l<(Ut/a)<t m p (7) 

where t m p is the time of maximum penetration. This time is specified in Eq. 8 below. 

TERMINAL REGIME 

The growth in inertial regime is arrested by strength and gravitational forces, and 
this marks the beginning of the terminal regime. In the case of simple craters (Figure 1) , 
the final shape of the crater is reached at the end of the inertial regime with the exception of 
the collapsing and folding over of the crater lip and minor elastic rebounding. In the case of 
complex craters (Figure 2), the end of the inertial regime marks the beginning of a series of 
complex motions which subsequently occur. These motions include the rebounding of the 
crater floor and the formation of central peaks and the collapse and propagation of the crater 
lip to form ringed craters. 
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The beginning of the terminal regime for crater depth occurs at the time of 
maximum penetration, tmp. This time, in the case where the gravitational forces dominant 
over the strength forces, is given by 

Ut,np/ a = 0.92 ( ga/U 2 (8) 

where the exponent magnitude is given by (l+p)/(2+p). The proportionality constant is 
obtained by fitting the results shown in Fig. 3. 

The maximum in the depth of penetration in the gravity dominated case is a 
function of the inverse Froude number and is obtained from Figure 5. The associated 
scaling law is given by 

Dp/a - 1.2 (ga/U 2 )-°- 22 (9) 

where the exponent is given by p/(2+p) and the proportionality contant is obtained from 
the results shown in Figure 5. 

The beginning of the terminal regime or the time of maximum penetration tmp, in 
the case where strength forces are dominant is given by 

Ut mp /a = 0.33 ( Yeff/pU 2 )-°-78 (10) 

where the exponent is (1 +p)/2 Here Y e ff is equal to twice the maximum shear strength. 
For the crusts of typical cratered planetary surfaces, we expect Y e ff will be significandy 
diminished from data from tests on undamaged rock on account of shock-induced cracking 
(e.g. Simmons et al. [1973] and Ahrens and Rubin [1992]). The temporal evolution in the 
strength dominated regime is shown in Fig. 6 from which Eq. 10 is derived by fitting the 
peaks of D/a versus Ut/a relations shown in Fig. 6. 

The maximum depth of penetration in the strength dominated case is a function of 
the inverse Cauchy number and is shown in Figure 7. The depth follows a power law 
with an expected slope given by p/2 = 0.28 until the inverse Cauchy number approaches ~ 
10 -2 . This change in scaling is a result of increased importance of shock weakening of the 
planetary material in the cratering process. Weakening is a result of the irreversible work 
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done in by the shock and subsequent unloading process and as a result of the post-shock 
deformation of the rock surrounding the crater. The irreversible work will result from both 
thermal shock heating and mechanical fracturing effects and is more intense close to the 
impact point where the high shock stresses are experienced. The importance of shock 
weakening depends upon the planetary strength, thermal state, and impact velocity. In the 
case of weak planetary crusts, and low impact velocities, the volume of the final crater is 
much greater than the shock weakened region and the effect is minor. However in the case 
of very strong initial crustal strengths, high surface temperatures or high (>30 km/sec) 
impact velocities, the extent of the shock weakening is comparable to the crater size, and 
the shock weakening effect can dominate late crater growth and the entire size and shape of 
the final crater. 

Shock weakening was modeled as thermal weakening and did not account for 
fracturing, which is currently an area of active interest (e.g.[Ahrens and Rubin, 1992; 
Asphaug et al., 1991; Housen and Holsapple, 1990]) and the present work should be 
extended to include that effect Our code computed the temperature field at each time step. 
The local temperature was used to compute the degree of thermal weakening by reducing 
the local strength in proportion to the difference between temperature and the melt 
temperature where the strength vanishes. The results of the code runs are plotted in Fig. 7 
for both depth of penetration and diameter. 

In the strength dominated case where shock weakening is negligible, the scaling is 
given by 

Dp/a = 0.84 ( Y/pU 2 J- 0 - 28 (11) 

which is derived from Fig. 7. 

In the shock weakened case, the depth is independent of the magnitude of the yield 
strength and is determined by the amount of material that is thermally weakened. In this 
case, the depth of penetration scaling is similar to depth of melting scaling [Bjorkman, 

1983] and is given by 
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Dp/a = 0.70 ( Hm/U 2 J* 0 - 28 ( 1 2) 

where Hm is the enthalpy required to induce shock melting, which in the case of thermal 
weakening is equal to the enthalpy to melt as referenced to the planets’ ambient surface 
condition. The constants in Eq. 12 are obtained from fitting the calculational results of Fig. 
8. The exponent is p/2, which is identical to the strength case. 

The condition for the transition between the strength and the shock weakened 
regimes can be obtained by equating Eqs. 1 1 and 12, which gives 

Y/pH m = 0.25 (13) 

A general relationship for depth that accounts for both regimes can be obtained by 
adding Eqs. 1 1 and 12 which yields 

Dp/a = K s ( Y/pU 2 )*0.28 +K SW ( Hm/U 2 (14) 

where K s = 0.84 and K sw = 0.7. This relationship can be plotted over a range of values 
relevant to planetary surfaces and is shown in Fig. 8. Eq. 14 implies an effective strength 
which accounts for shock weakening which can be defined as 

Y ef f = Y { 1 + (K sw /K s ) (Y/pH m )0.28}(- 1/0.28) (15) 

This result is plotted in Figure 9. 

A general interpolation relationship of the form suggested by Holsapple and 
Schmidt [1987] can be used to span the gravity and strength regimes 

Dp/a = K s (ga/U 2 )-0- 22 { (K g /K s )-0- 22 + (Y/pU 2 )^ 22 ^ 28 ) ( ga /lJ2)-i }-0.22 (16) 
where Kg = 1.2. 

Replacing the shear strength by the effective strength defined in Eq. 15 
leads to a more general scaling relationship that includes gravity, strength and shock 
weakening effects 

Dp/a = K s (ga/U 2 )-°- 22 {(K g /K s )-0- 22 + { (Y/pU 2 ) { 1 + 

(K sw /K s ) (Y/pH m )028 } .1/6.28) (0.22/0.28 ) (ga/U 2 )-!}^ 22 } (17) 

The maximum depth of penetration of projectiles from giant impacts has been a 
topic of much debate in the literature (e.g. Grieve et al.[1981]). The tracer particles 
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employed in the calculations help elucidate the motions of material in the cratering process. 
Referring to Figs. 1 and 2, we see that near the centerline of impact, planetary material is 
driven downward and is not ejected in the cratering process. As one moves from the 
centerline towards the crater lip, there is a region in which the planetary material is 
essentially stagnated against the outwardly moving cavity wall; this stagnated region, 
delineated by layers of tracer particles, shows a lower layer of particles being driven 
downward and an upper layer moving upward and being excavated (see Fig. lc). The 
scaling laws for the maximum depth of excavation is given by 

Dex/a = 0.17 (Y/pU 2 )* 0 - 28 (18a) 

and 

Dex/a=K(ga/U2)-0- 22 (18b) 

for the strength dominated regime. Here the exponents for strength and gravity scaling are 
given by Eqs. 1 1 and 12. The proportional constant is obtained from Figs. 2 and 3. 

Pi a m ft CT 

The beginning of the terminal regime for the diameter occurs when the crater lip 
stops expanding laterally. In the case where the gravitational forces are dominant, this 
occurs when 

Ut/a = 1.8 ( ga/U 2 )’®- 61 (19) 

where the exponent magnitude is given by (l+p)/(2+|i) and the proportionality factor in 
this equation and those of Eqs. 20 and 21, below, come from fitting the results of Fig. 4. 

An example of the crater morphology at the time the crater cavity has reached its 
maximum diameter in the gravity dominated case is shown in Figure 2c. The cavity 
maximum diameter is a function of the inverse Froude number and is given by 

dm/a = 2.1 (ga/U2)-0- 22 (20) 

where the exponent is p/(2+p). 

At the time the crater cavity has reached its maximum diameter, the crater lip has 
reached its maximum height (figure 2c). This height scales as 
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H*/a= 0.53 ( ga/U 2 ) -0 - 22 


( 21 ) 


where the exponent is again p/(2+p)where the proportionality factor was obtained from 
plots such as Fig. 2c. 

The crater lip will collapse and produce an outwardly propagating surface wave 
which rapidly increases the apparent diameter of the crater (Figs. 2c-2d). Simultaneously, 
the crater floor moves upward and an inner crater lip forms. Note that neither the outer ring 
nor the inner ring correspond to the maximum diameter of the transient crater. 

The diameter of the outer ring scales as 

d/a = 4.28 ( ga/U 2 )- 0 - 22 (22) 

where as before the exponent is given by |i/(2+p). The proportionality factor comes from 
Eq. 4. 

When strength dominates the cratering process, the crater growth ends in the inertial 
regime with only a small amount of elastic rebounding occurring followed by long term 
crater relaxation. The time of termination of crater formation is given by 

Ut/a = 0.33 ( Yeff/pU2)-0.78 (23) 

Note that this time is identical to that for the depth. The parameters come from fitting the 
results of Fig. 6. 

The maximum or final crater diameter in the strength dominated regime is given by 

d/a =1.4 (Y/pU 2 )- 0 - 28 (24) 

Here we employed the results shown in Fig. 7. In the shock weakened regime the final 
diameter is given by 

d/a= 1.3 (Hm/U 2 )- 0 - 28 (25) 

where the exponent is p/2 . 

In a similar manner to the depth, a general interpolation relationship spanning the 
strength and shock weakened regimes is given by 

d/a = K s (Y/pU2)-0-28 + Ksw (Hm/U 2 )" 0,28 , (26) 
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where for the diameter K s = 1.4 (Eq. 24), and K sw =1.3 (Eq. 25). 


Transition from Simple to Complex Craters 

The transition from simple to complex craters is the transition from strength- 
dominated to gravity-dominated craters. A general interpolation relationship for the crater 
diameter which spans both the gravity and strength regimes is given by 
d/a = K s (ga/U2)-0 22 { (K g /K s )-°-22 + ( (Y/pU 2 ){l + 

(K sw /K s ) (Y/pH m ) 0 - 28 )-lA)-28} (0.22/0.28) (ga AJ 2 )' 1 }’ 0 22 (27) 

where for the diameter K s = 1.4 (Eq. 24), and K sw =1.3 (Eq. 25), and Kg = 4.28 (Eq. 
22). This result is plotted in Figure 10. 

The transition diameter, or equivalently, the intersection between the strength and 
gravity dominated regimes can be readily derived . The transition diameter is obtained by 
equating Eqs. 22 and 24, and solving for the transition impactor radius. Substituting this 
radius into Eq. 22 yields the following simple relationship 

d t = (K d g /Kd s )((2+p)/W) Y/pg (28) 

where K d g is the gravity proportionality constant in Equation 23 and K d g is the strength 
proportionality constant in Equation 25, and using those numerical values we obtain 
(K d g/K d s )((2+ti)/ti)) = 8.6. Note that the transition diameter is directly proportional to the 
strength and inversely proportional to gravitational acceleration. This result is compared in 
Fig. 1 1 to the transition diameters as a function of surface gravity for the terrestrial planets 
and various satellites [Chapman and McKinnon, 1986; Schenk, 1989]. 

The transition from simple to complex craters has been studied for some time using 
static models to predict the slumping of freshly formed craters [Melosh, 1977; Melosh and 
McKinnon, 1978]. Melosh [1977] found, using a perfectly plastic constitutive relation, 
that the criterion for stability of freshly formed bowl-shaped craters is given by 

pgDp/Y = K m 5< K m <10 (29) 
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We will show that Eq. 29 can be cast in the form of Eq. 28 and leads to similar 
conclusions as those drawn by Melosh. Eq. 29 can be solved for the maximum depth of 
penetration and that result is used in Eq. 1 1 to determine the transition impactor radius. 
Substituting this into Eq.24 for the diameter yields 

dt= (K m K^s/K D s ) Y/pg (30) 

where K D S is the strength proportionality constant in Eq. 11. Equating Eq. 28 and 30 gives 

the following expression for the Melosh constant 

K m = KD s (Kd g d l /Kd s )((2+p)/u) (31) 

For the numerical values found in our calculational results, we obtain a value of K m = 8.9, 
which is in agreement with Melosh’s range of values (5=>10). 

Shape . Observables, and Inferables 

In studying cratering on remote planetary surfaces, the observables are the surface 
morphology of the craters. In the following section we will relate some of the transient 
non-observable aspects of the cratering process (e.g. depth of penetration, excavation ....) 
to the usual observable quantity — crater diameter. 

The maximum transient cavity diameter relative to the final crater diameter in gravity 
dominated craters is given by dividing Eq. 20 by Eq. 22, which yields 

d„»/d = 0.5 (32) 

and for strength dominated craters, the maximum cavity diameter is equal to the final 
diameter, except for minor crater rebounding. Notably, the maximum transient diameter is 
the maximum diameter in this case, 

dtm/d = 1 (33) 

The maximum crater lip height relative to the crater diameter in gravity dominated craters 
can be obtained by dividing Eq. 21 by Eq. 22 to get 

H//d = 0.12 (34) 
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The ratio of the time for the crater lip to collapse (tip) to the time for the crater floor 
to rebound (tf) in the gravity dominated regime can be obtained from Figure 4 and is given 
by 

t|p/tr = 2 (35) 

The depth of penetration relative to the crater diameter for gravity-dominated 
craters can be obtained by dividing Eq. 9 by Eq. 22, which yields 

Dp/d = 0.28 (36) 

and for strength dominated craters by dividing Eq. 1 1 by Eq. 24, which yields 

Dp/d = 0.58 (37) 

In a similar manner, the depth of excavation relative to the crater diameter for 
gravity dominated craters can be obtained by dividing Eq. 18a by Eq. 22, which yields 

D ex /d = 0.1 (38) 

and for strength dominated craters by dividing Eq. 18b by Eq. 24 

D ex /d = 0.2 (39) 

Note that the depth of penetration and excavation relative to the crater diameter are a 
factor of two less in the gravity dominated cases because the crater diameter is larger as a 
result of the lip collapse and lipwave propagation. The above result for the depth of 
excavation is compared to field data from terrestrial craters in Fig. 12. The maximum depth 
of penetration and depth of excavation as a function of crater diameter are plotted for the 
various planets in Fig. 13. As an example, for craters in the range of 150 km in diameter, 
which is in the estimated size for the K-T event [Alvarez et al., 1980], the maximum depth 
of penetration is 42 km while the maximum depth of excavation is -15 km. Impacts of 
this magnitude appear to open the possibility of impact-induced volcanism from the mantle. 


CONCLUSIONS 

The objective of this study was to obtain a quantitative understanding of the 
cratering process over a broad range of conditions. In accomplishing this objective we 
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numerically computed the evolution of impact induced flow fields and calculated the time 
histories of the major measures of crater geometry (e.g. depth, diameter, and lip height) for 
variations in planetary gravity, material strength, thermal properties and impactor size. The 
results of these calculations are summarized in Table 2. These results were fit with the 
scaling laws of Holsapple and Schmidt [1987]. In all of the measures we examined, these 
scaling laws were found to be in good agreement with the numerical results. From this 
experience, we expect this agreement to hold in general. 

We have described the impact process in terms of four regimes: (1) penetration, (2) 
inertial, (3) terminal , and (4) relaxation. The scaling laws for each of the first three 
regimes are given in Table 3. In the penetration regime the depth evolves linearly with 
time. In the inertial regime, the scaling of the crater geometry is a function of the coupling 
coefficient (|i) as defined by Holsapple and Schmidt [1987] and was shown by them to be 
related to the slope of the growth of the depth and diameter in that regime. The results of 
this paper confirm that relationship. 

In studying cratering on remote planetary surfaces, the primary observables are the 
crater diameter including the diameters of the multiple rings, the presence of central peaks 
and pits and the relaxed crater depth. We have related several key transient non-observable 
aspects of the cratering process (e.g. maximum depth of penetration, depth of excavation, 
and maximum crater lip height) to the crater diameter. For example, the maximum depth of 
penetration relative to the crater diameter is ~0.58 for strength dominated craters and -0.28 
for gravity dominated craters. This implies that impactors associated with the large basin 
impacts penetrated deep into the planetary surfaces. In examining the details of the 
rebounding of the depth in the gravity dominated regime, we found that the initial planetaiy 
material layering sequence is reconstructed in the rebound process. This phenomenon thus 
accounts for apparent contradictions of field data and experiments conducted on the 
geotechnical centrifuge [Schmidt and Holsapple, 1981], that the impactor did not penetrate 
very deep because the layering sequence was relatively undisturbed. In a similar manner 
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to the penetration, the maximum deep of excavation relative to the crater diameter is ~0.1 
for gravity dominated craters and -0.2 for strength dominated; this implies that for impact 
velocities less than the threshold for significant vaporization, the excavated material can 
have rather deep origins for large diameter basin-forming craters 

The transition from simple bowl shaped craters to complex craters was analyzed. 
This occurs at the transition from strength to gravity dominated craters. The transition 
diameter is proportional to Y/pg where the proportionality constant is -8.6. This constant 
is in the range of values (5-10) found by Melosh [1977] using a gravitational plastic 
instability model. 
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Table 1. Scope of Parameters Studied for Impact of 
Silicate Projectile on a Planetary Halfspace. 


Parameter 

Symbol 

Values Employed 

Units 

impact velocity 

U 

12 

km/sec 

planetary gravity 

g 

O.UO^IOMO 6 

ge (980cm/s 2 ) 

density 

P 

2.7 

g/cm^ 

bulk modulus 

E 

7.6 x 10 11 

dynes/cm 2 

Gruneisen coefficient 

Y 

2.0 

— 

melt enthalpy 

H m 

1.1 x 1010 

ergs/g 

vaporization enthalpy 

Hvap 

1.2 x 101 1, 1.2 x 10*2, 
8.0 x 10 12 

ergs/g 

yield strength 

Y 

0,0.24,2.4,5.6, 24, 28,140, 240,2400 

kbar 

inverse Froude Number 

ga/U 2 

0,3.4x10-7, 3.4x10*5, 3.4x10-3, 0.34 

— 

inverse Cauchy Number 

Y/pU 2 

0, 7.2x10-5,7.2x10-^,7.2x10-3,0.036 

— 

melt number 

Hm/U 2 

7.6 x 10-3 


vaporization number 

Hvaj/U 2 

0.08, 0.8, 1.2 






Table 2. Normalized Crater Diameter and Depth, versus, x 
for different values of ga/U 2 , and Y/pU 2 . The impact velocity, U, 
was 12 km/s except where indicated. 


ga/U 2 Y/pU 2 


0 0.0 

" 0.0 

" 0.0 

0 (U=10km/s) a 0.0 

" 0.0 

" 0.0 

M 0.0 

" 0.0 

0 (U=50km/s) a 0.0 

" 0.0 

" 0.0 

” 0.0 


" 0.0 

" 0.0 

0 (U=100km/s) a 0.0 

" 0.0 

" 0.0 

" 0.0 

’’ 0.0 

" 0.0 

" 0.0 

” 0.0 

0.0 

" 0.0 

3.41 xlO- 7 1.4x10*3 

1.4x10*3 
1.4x10*3 
1.4x10*3 
1.4x10*3 

" 1.4x10*3 

1.4x10*3 
1.4x10*3 
1.4x10*3 
1.4x10*3 
1.4x10*3 

” 1.4x10*3 

1.4x10*3 

3.41xl0- 7 7.2x10*3 

7.2x10*3 


d/a 

D/a 

X 


0.43 

0.870 


1.07 

2.140 


1.53 

3.119 


0.2 

0.4 


0.8 

1.6 


2.0 

5.0 


3.0 

12.0 


4.0 

16.0 


0.4 

0.8 


1.0 

2.0 


2.0 

4.0 


2.6 

6.0 


4.0 

19.0 


5.2 

50.0 


0.2 

0.41 


1.0 

2.00 


2.0 

4.00 


2.6 

6.00 


4.0 

20.00 


5.0 

40.00 


6.0 

75.00 


7.0 

105.00 


8.0 

130.00 


9.5 

200.00 

3.400 

1.65 

3.18 

5.100 

3.00 

9.28 

5.700 

3.35 

12.28 

6.500 

3.75 

16.50 

6.500 

3.75 

16.50 

6.600 

3.95 

18.50 

7.400 

4.40 

23.30 

7.800 

4.75 

29.64 

8.800 

5.15 

41.00 

9.000 

5.20 

55.90 

9.000 

5.10 

75.40 

9.000 

5.10 

94.00 

9.000 

5.00 

113.00 

3.500 

1.700 

3.480 

4.500 

2.700 

7.700 



Table 2. Normalized Crater Diameter and Depth, versus, X 

for different values of ga/U 2 , and Y/pU 2 . The impact velocity, 
U, was 12 km/s except where indicated.(continued) 


ga/U 2 

Y/pU2 

d/a 

D/a 

X 

*i 

7.2x10*3 

5.000 

3.000 

10.560 

•t 

7.2x10-3 

5.500 

3.500 

14.800 

tt 

7.2x10*3 

5.700 - 

3.700 

17.600 

it 

7.2x10*3 

6.000 

3.750 

22.100 

•i 

7.2x10*3 

6.000 

3.600 

26.400 

•t 

7.2x10*3 

5.700 

3.400 

30.480 

•t 

7.2x10*3 

5.600 

3.200 

34.300 

tt 

7.2x10*3 

5.400 

2.800 

37.400 

3.41xl0- 7 

3.6xl0- 2 

3.300 

1.750 

3.360 

tt 

3.6x10*2 

4.350 

2.900 

8.900 

tt 

3.6x10*2 

4.600 

3.100 

13.300 

tt 

3.6x10-2 

4.600 

2.950 

17.200 

tt 

3.6x10*2 

4.500 

2.700 

21.020 

ti 

3.6x10*2 

4.200 

2.500 

24.400 

tt 

3.6x10-2 

4.000 

2.450 

26.800 

tt 

3.6x10*2 

3.800 

2.300 

29.500 

tt 

3.6x10*2 

2.100 

2.200 

34.080 

tt 

3.6x10-2 

1.000 

1.500 

49.600 

3.41x10-5 

0.0 


0.37 

0.825 

tt 

0.0 


1.16 

2.350 

tt 

0.0 


2.17 

5.590 

tt 

0.0 


2.60 

7.700 

tt 

0.0 


1.67 

3.600 

tt 

0.0 


3.02 

10.100 

tt 

0.0 


3.53 

14.180 

it 

0.0 


4.05 

19.000 

tt 

0.0 


4.42 

24.000 

3.41xl0- 3 

0.0 


0.40 

0.825 

tt 

0.0 


1.90 

4.480 

tt 

0.0 


2.63 

7.890 

tt 

0.0 


3.25 

13.279 

tt 

0.0 


4.10 

22.770 

tt 

0.0 


4.20 

32.000 

tt 

0.0 


3.10 

41.200 

tt 

0.0 


1.40 

50.600 

it 

0.0 


1.20 

2.360 

tt 

0.0 


4.10 

27.600 

it 

0.0 


3.45 

36.800 

3.41x10-3 

6.2x10-5 


0.75 

1.550 

tt 

6.2x10-5 


0.40 

0.854 

tt 

6.2x10-5 


1.20 

2.370 
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Table 2. Normalized Crater Diameter and Depth, versus, x 
for different values of ga/U 2 , and Y/pU 2 . The impact velocity, 
U, was 12 km/s except where indicated.(condnued) 


ga/U 2 

Y/pU 2 

d/a 

D/a 

X 

tt 

6.2x10*5 


1.50 

3.090 

it 

6.2x10-5 


1.60 

3.370 

« 

6.2x10*5 


1.70 

3.640 

it 

6.2x10*5 


2.00 

5.200 

it 

6.2x10*5 


2.30 

6.070 

tt 

6.2x10-5 


0.10 

0.200 

ft 

6.2x10-5 


2.70 

8.350 

tt 

6.2x10-5 


2.50 

7.200 

tt 

6.2x10-5 


2.90 

9.660 

tt 

6.2x10-5 


3.20 

12.690 

tt 

6.2x10-5 


3.60 

17.100 

tt 

6.2x10*5 


3.80 

21.900 

tt 

6.2x10-5 


3.90 

27.800 

tt 

6.2x10-5 


3.65 

35.400 

tt 

6.2x10-5 


3.00 

42.600 

ti 

6.2x10-5 


2.00 

50.400 

3.41x10-3 

6.2x1 0- 4 


0.50 

0.907 

tt 

6.2x1 0 -4 


1.10 

2.052 

tt 

6.2x1 0- 4 


1.70 

3.280 

tt 

6.2x1 O* 4 


2.15 

4.776 

3.41x10*3 

6.2x10-3 


0.50 

0.934 

it 

6.2x10*3 


1.22 

2.324 

tt 

6.2x10-3 


1.85 

3.710 

tt 

6.2x10*3 


2.35 

5.560 

tt 

6.2x10-3 


2.55 

7.040 

tt 

6.2x10*3 


3.00 

10.150 

tt 

6.2x10-3 


3.35 

15.090 

tt 

6.2x10-3 


3.30 

20.020 

tt 

6.2x10*3 


2.20 

43.200 

tt 

6.2x10*3 


1.85 

50.850 

tt 

6.2x10-3 


3.00 

27.600 

tt 

6.2x10-3 


2.50 

35.470 

3.4xl0- 2 

0 

1.250 

1.400 

21.670 

tt 

0 

0.900 

0.600 

28.800 

tt 

0 

0.700 

0.200 

33.600 

tt 

0 

0.250 

0.100 

37.700 

tt 

0 

2.450 

0.100 

42.960 

tt 

0 

4.100 

0.100 

46.560 

3.4xl0- 2 

6.2x10-5 

2.000 

0.400 

0.828 

tt 

6.2x10-5 

2.750 

1.100 

2.120 


3 
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Table 2. Normalized Crater Diameter and Depth, versus, t 
for different values of ga/U 2 , and Y/pU 2 . The impact velocity, 
U, was 12 km/s except where indicated.(continued) 


ga/U 2 

Y/pU 2 

d/a 

D/a 

t 

tt 

6.2x10-5 

3.300 



H 

6.2x10-5 

3.550 



»l 

6.2x10-5 

3.850 



fl 

6.2x10*5 

4.100 



M 

6.2x10*5 

4.250 



•1 

6.2x10*5 

4.850 



It 

6.2x10*5 

6.600 



It 

6.2x10-5 

8.350 



•1 

6.2x10*5 

8.600 



It 

6.2x10*5 

8.100 



•I 

6.2x10*5 

8.200 



•1 

6.2x10*5 

8.600 



tl 

6.2x10*5 

8.600 



3.41xl0- 2 

6.2X10- 4 

5.200 



it 

6.2x10-* 

9.150 



it 

6.2xl0 4 

9.100 



it 

6.2x10-4 

9.100 



tt 

6.2x10-* 

9.100 



it 

6.2x10-4 

9.100 



tt 

6.2x10-* 

10.000 



it 

6.2x10-4 



49.920 

3.4xl0- 2 

6.2x10*3 

2.700 

1.350 


tt 

6.2x10-3 

2.950 



it 

6.2x10-3 

4.400 



tt 

6.2x10-3 

4.600 



tt 

6.2x10*3 

5.050 



it 

6.2x10-3 

5.550 



it 

6.2x10-3 

6.500 



•t 

6.2x10*3 

7.500 



it 

6.2x10*3 

8.200 



•t 

6.2x10*3 

9.400 



3.4x10*2 

6.2x10*2 

4.850 



it 

6.2x10-3 

5.100 



it 

6.2x10*3 

5.300 



it 

6.2x10*3 

7.500 



tt 

6.2x10-3 

8.100 



3.4x10-2 

6.2x10-1 

3.000 



ft 

6.2x10-1 

3.350 



it 

6.2x10-1 

3.800 
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Table 2. Normalized Crater Diameter and Depth, versus, x 

for different values of ga/U 2 , and Y/pU 2 . The impact velocity, 
U, was 12 km/s except where indicated.(condnued) 


ga/U 2 

Y/pU 2 

d/a 

D/a 

T 

ft 

6.2x10-1 

4.000 

2.400 

8.640 

ft 

6.2x10-1 

3.950 

2.350 

9.980 

it 

6.2x10-1 

3.950 

2.350 

10.500 

it 

6.2x10-1 

4.150 

2.200 

12.600 

it 

6.2x10*1 

4.200 

1.950 

18.000 

3.4xl0- 2 

6.2 

3.200 

1.450 

2.660 

it 

6.2 

3.800 

2.400 

6.766 

•f 

6.2 

4.050 

2.400 

9.040 

it 

6.2 

4.150 

2.350 

11.400 

it 

6.2 

4.200 

2.250 

12.900 

it 

6.2 

4.050 

2.000 

16.300 

0.341 

0.0 

2.1 

0.40 

0.765 

it 

0.0 

2.8 

1.10 

2.010 

it 

0.0 

2.3 

1.50 

4.080 

tt 

0.0 

5.0 

1.35 

5.220 

tt 

0.0 

5.8 

1.10 

5.970 

it 

0.0 

2.6 

1.50 

3.500 

it 

0.0 

2.0 

1.40 

4.480 

it 

0.0 

6.8 

1.10 

7.050 

it 

0.0 

2.4 

1.50 

3.900 

it 

0.0 

2.3 

1.45 

4.180 


0.0 

4.0 

1.30 

4.750 


0.0 

5.3 

1.25 

5.678 

0.341 

6.2x10*5 


0.38 

0.761 

•t 

6.2x10*5 


0.68 

1.290 

it 

6.2x10-5 


0.90 

1.700 

tt 

6.2x10-5 


1.15 

2.130 

•t 

6.2x10-5 


1.30 

2.490 

it 

6.2x10*5 


1.35 

2.950 

it 

6.2x10*5 


1.43 

3.400 

it 

6.2x10*5 


1.50 

2.880 


(a) Time history taken from Holsapple & Schmidt (1987). 
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TABLE 3. SUMMARY OF CRATER SCALING FORMULA 


Regime 

Formula 

Proportionality 
Constant (K) 

Exponent 

Penetration Regime 

Depth Evolution 

D/a = 0.5(Ut/a) 

0.5 

1 

Depth transition to inertial 

Ut/a=5.1 

— 


regime 

Diameter Evolution 
Diameter transition to inertial 

Ut/a=2 

1.3 

0.78 

regime 

Diameter/depth evolution 

d/D=K(Ut/a) 



Inertial Regime 

Depth Evolution 

D/a=K(Ut/a) s 

1.3 

0.36 

Coupling coefficient 

\i = s/(l-s) 

0.56 

— 


Transition to terminal (g)* 

Ut/a=K(gaAJ2)( 1+ PX2+p) 

0.92 

0.61 

Transition to terminal (Y) + 

Ut/a=K(Y/pU 2 )-( 1+ P>' 2 

0.33 

0.78 

Maximum depth of penetration (g) 

Dp/a = K (ga/U 2 )-^ 2+ P) 

1.2 

0.22 

Maximum depth of penetration (Y) 

Dp/a=K(Y/pU 2 )-M^ 

0.84 

0.28 

Maximum depth penetration(Y sw ) ++ 

D p /a=K(H m AJ 2 )-^ 

0.70 

0.28 


Comments 


Independent of 
material properties 


l<(Ut/a)<5.1 


Independent of 
material properties 
Slope of growth, s, 
determines the value 
of (i = s/(l-s) 
which, in turn, 
determines the 
exponents for all 
regimes (Holsapple 
and Schmidt, 1987) 


TABLE 3. SUMMARY OF CRATER SCALING FORMULA (continued) 


Penetration Regime 

Formula Constant (K) 

Exponent 

Comments 

Diameter Evolution 

d/a=K (Ut/a) s 

2.0 

0.36 

Crater grows until 

Maximum cavity diameter (g) 

d m /a=K(ga/U 2 )-^ 2+ P) 

2.1 

0.22 

arrested by strength 

Maximum diameter (Y) 

d m /a=K(Y/pU 2 )-^ 

1.4 

0.28 

and/or gravitational 

Maximum diameter (Y sw) 

d m /a=K(H m /U 2 r ,jy2 

1.3 

0.28 

forces 

Transition time to terminal (g) 

Ut/a=K(ga/U 2 )*< 1+ M 2+ P) 

1.80 

0.61 


Transition time to terminal (Y) 

U t/a=K( Y/pU 2 )"^ 1 + PX 2+ H) 

0.33 

0.78 


Maximum lip height (g) 

H//a=K(ga/U 2 )*M/( 2+ P> 

0.53 

0.22 


Piameter/Depth Evolution 

d/I^KOlt/aJ-^O+P) 

d/D=K 

4.0 

1.6 

0.64 

l<(Ut/a)<5.1 
5. 1 <(Ut/a<)tfnp 

Terminal Regime 

Transition diameter from 
simple to complex craters 

d t =K(Y/pg) 

8.6 

— 

In strength 
dominated regime 

Maximum diameter (g) 

Maximum transient diameter/ 
diameter (g) 

Maximum transient diameter/ 
diameter (Y) 

Maximum lip height/diameter 
Depth of penetration/diameter (g) 
Depth of penetration/diameter (Y) 
Depth of excavation/diameter (g) 
Depth of excavation/diameter (Y) 

dm/a = K ( ga /U2)-M/(2+4) 

dm/d — 0.5 

dm/d = 1.0 
H/ d = 0.12 
Dp/d = 0.28 
Dp/d = 0.58 
D ex /d = 0.1 
Dex/d = 0.2 

4.28 

0.22 

only minor elastic 
rebounding . In 
gravity dominated 
regime crater 
evolves until 
arrested by strength 
forces 


* (g) Regime indicated is gravity dominated 

+ (Y) Regime indicated is strength dominated 

++ Ysw Regime indicated is shock wave degraded strength dominated 


2 




Figure 1 . Simple bowl shaped crater Flow field for the impact of a silicate projectile on a 
silicate halfspace. The flow field is in the strength dominated regime with ga/U 2 =3.7x10*7 
and Y/pU^ 1.4x1 0*3. The velocity field is shown on the left side with the scale shown at 
the top. On the right hand are the trajectories of tracer particles placed at various depths. 
The panels represent dimensionless times of 3.19, 9.29 and 16.6. This is representative of 
a simple bowl shaped crater. 

Figure 2. Complex crater Velocity and deformation history fields for the impact of a 
silicate projectile on a silicate halfspace, the flow field is gravity dominated at late times 
with ga/LK = 3.4 x 10* 2 and Y/pU^ = 6.2 x 10*3. The velocity field is shown on the left 
side with the scale at the top. On the right side are the trajectories of tracer particles placed 
at various depths. The panels represent dimensionless times of 5.57, 8.32, 10.2, 12.8, 
14.4, 17.9, and 19.7. 

Figure 3. Dimensionless depth versus dimensionless time for the penetration, and inertial, 
terminal and relaxation regimes for different values of ga/U 2 and Y/pU 2 

Figure 4. Temporal evolution of the dimensionless depth of penetration (D/a) and crater 
diameter (d/a) for a fixed value of inverse Froude number (ga/U 2 =0.034), and various 
values of Y/pga. 

Figure 5. Dimensionless maximum depth of penetration (Dp/a) and diameter (d/a) as a 
function inverse Froude number (ga/U 2 ). Slope corresponds to - p/(2+p) for gravity 
scaling. Symbols are results for different calculations. 

Figure 6. Temporal evolution of dimensionless depth (D/a) and diameter (d/a) in the 
strength dominated region. 

Figure 7. Dimensionless maximum depth of penetration (Dp/a) and diameter (d/a) as a 
function of inverse Cauchy number (Y/pU 2 ). The slope is -0.28=-|i/2 and corresponds to 
strength scaling. The change in slope is a result of melting and shock weakening. 

Symbols represent different calculations for a fixed value of Froude number (ga/U 2 =0.034) 
for various values of the melt number (Hm/U 2 ). 

Figure 8. Dimensionless crater diameter (d/a) as a function of inverse Cauchy number 
(Y/pU 2 ) in the strength dominated regime (see Figure 10). 

Figure 9. Effective strength as a function of weakening energy. 

Figure 10. Dimensionless crater diameter as a function of inverse Froude number (ga/U 2 ) 
for various values of inverse Cauchy number (Y/pU 2 ). Slope corresponds to gravity 
scaling. 

Figure 1 1. Crater transition diameter from simple to complex craters as a function of 
planetary surface gravity. Slope corresponds to scaling law of Eq. 28. 

Figure 12. Depth of excavation of planetary material (D ex ) divided by the Crater diameter 
(d) as a function of crater diameter for terrestrial craters. Calculation results indicated by 
solid line, (a) sand, Gault et al., 1968; (b) sand, Andrews, 1976; (c) 500 ton TNT, 
Roddy, 1976; (d) Jangle U, Shoemaker, 1963; (e) Teapot Ess, Shoemaker, 1963; (f) 



Odessa, Shoemaker and Eggleston, 1961; (g) LonarLake, Fredriksson et al., 1973; (h) 
Decaturville, Offield and Pohn, 1976. 

Figure 13. Maximum depth of penetration (Dp) and excavation (D e x) as a function of crater 
diameter for various planets and satellites. Vertical lines delineate the transition from 
strength to gravity scaling. 
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